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I. INTRODUCTION 



Recently, the production of high-density and high-temperature matter in relativistic nu- 
clear collisions has received much attention. It is generally expected that these collisions 
may provide the tools to probe the existence of a new phase of matter of strongly interacting 
particles, the quark-gluon plasma (QGP), at high temperature or at high baryon density 
[1]. In the QGP, quarks and gluons are deconfined; the deconfinement refers to a circum- 
stance in which quarks, anti-quarks, and gluons are no longer confined within the spatial 
dimensions of a hadron. Quarks and gluons are nonetheless confined within the boundary 
of the plasma. Under laboratory conditions in which a QGP is produced using relativistic 
heavy-ion collisions, the expected initial transverse radius of the QGP is comparable to the 
smaller of the radii of the two colliding nuclei producing the plasma, typically a few fermis. 
Therefore, a corollary of the existence of a QGP in a heavy-ion collision is that the quarks 
and gluons in the plasma travel freely within a spatial region of a few fermis. 

The distribution of transverse momentum provides information on the transverse bound- 
ary of the plasma. It can be used to detect the presence of the QGP. Furthermore, the 
magnitudes of the signals for QGP detection are functions of the momentum distribution of 
the quarks and gluons, which has generally been taken to be a thermal distribution [2-4]. 
Because the presence of the boundary affects the momentum distribution [5], the boundary 
effects should be taken into account to give better estimates of the magnitude of the sig- 
nals of the quark-gluon plasma. While the boundary effects on quarks have been studied 
previously [5], here we wish to study the effect of the boundary on the distribution of the 
transverse momenta of gluons in a QGP, using wave functions which satisfy the boundary 
conditions. 

Another interesting boson system in which boundary effects may manifest themselves 
is the multipion system produced after a high energy heavy-ion collision. If the mean-field 
potential of the pion system has a boundary before it freezes-out, the transverse momentum 
distribution will show the effect of the boundary. 

The subject of the transverse momentum distribution of pions has received much at- 
tention recently. The latest experimental data [6-8] show that the transverse momentum 
distributions of negative hadrons, mainly pions, from hadron-nucleus and nucleus-nucleus 
collisions exhibit a striking enhancement at low and high transverse momenta relative to 
that of negative hadrons from p -\- p collisions at comparable energies. This origin of this 
puzzling phenomenon has been vigorously debated in many articles [9-17]. Atwater et al. [9] 
and Lee and Heinz [10] interpreted the data by assuming a transversely expanding fireball. 
With assumed velocity profiles, Kusnezov and Bertsch [11] argued that the transverse veloc- 
ity profiles used in references [9,10] imply an unnatural assumption concerning the nature of 
the freeze-out surface. Furthermore, they indicated that a more complete hydrodynamical 
description of the transverse fiow cannot account for the experimental shape of the pion 
transverse momentum distribution at low p^. Another explanation, suggested by several 
authors [12-14], is a non-zero chemical potential for the pions. Alternatively, SoUfrank, 
Koch, and Heinz [15] found agreement with the S-|-S data by including contributions due to 
resonance decays on hadronic momentum spectra. However, Barz et al [16] examined the 
role of decays of the excited hadrons in determining the pion p^ spectra, and argued that 
no reasonable statistical model can reproduce the experimentally observed peak at low p^ 
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in the pion transverse momentum distributions. A very good agreement with the S+S data 
has also been found by Bolz et al [17] by modifying Landau initial conditions with an initial 
rapidity distribution which is only a function of the longitudinal coordinate. 

In this work we study the boundary and Coulomb effects on pion systems at high tem- 
perature. The momentum distribution of particles in a bounded medium is altered because 
the single-particle wave function must satisfy the appropriate boundary conditions. The 
pions produced after a nucleus-nucleus collision are moving in the Coulomb field of the col- 
liding nuclear participants, whose charge depends on the impact parameter. For a head-on 
collision, the charge of the colliding participants in the central rapidity region can be as big 
as 2 Z for two identical nuclei. The momentum distribution of charged pions is changed due 
to the Coulomb final-state interactions between the pions and the protons near the center 
of mass region. This Coulomb effect depends on the factor rj = Zajv, which is large in 
heavy-ion central collisions. This final state interaction has a large effect on the charged 
pions with small transverse momenta, as shown in this paper. 

This paper is organized as follows: In Section II, we describe the method for obtaining 
the momentum distribution of a boson system. In Section III, the transverse momentum 
distribution of gluons in a quark-gluon plasma is obtained, using wave functions satisfying 
the boundary conditions. In Section IV, we study the boundary effect on the pion transverse 
momentum distribution in a pion system at high temperatures. The effect of the Coulomb 
final-state interaction on the pion transverse momentum distribution is discussed in Section 
V. Finally, Section VI contains the summary and conclusions. 



For a boson system in a mean-field scalar potential m(x), if $a(p) is the boson single- 
particle wave function in a quantum state A in momentum space, the momentum distribution 
of bosons is given by: 



where the summation is carried out over all the single-particle quantum numbers {A}, 
is the boson degeneracy number, and W{E\) is the occupation probability of the single- 
particle state A. The thermodynamical state of the system is described by the distribution 
11^ (-Ea), which is a function of the single-particle energy E\. For a boson system at thermal 
equilibrium W{E) is represented by the Bose-Finstein distribution, which for a boson system 
at temperature T and chemical potential // is given by 



The wave function $(p), for a boson system with a cylindrical boundary, is given by the 
Fourier transform 



where $(r, ^,2;) is the corresponding single-particle wave function in configuration space, 
which is obtained by solving the Klein-Gordon equation for a boson in a mean-field potential 



II. MOMENTUM DISTRIBUTION OF A BOSON SYSTEM 





W/(^) = [e(^-^)/^-l]-i. 
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m(r). The Klein-Gordon equation, for the wave function $(r, ^,2;) of a boson in a boson 
system with a cyhndrical boundary and moving in a mean-field scalar potential m(r) varying 
in the transverse direction, is given by 

- m2(r)] $(r,(^,z) = 0, (2) 

where p is the boson four-momentum. Different boson systems will be described by different 
mean-field potentials m(r). 



III. BOUNDARY EFFECTS ON A GLUON PLASMA 



We consider a gluon in a quark-gluon plasma in the form of a sharp cylinder, which is 
presumed to be formed in a high energy nucleus-nucleus collision. The mean-field potential 
m(r) is defined by: 

, , j if r < Ro (region I) 
] M if r > Rn (region II) , 



mi r 



for gluons we set = 0. Separation of the wave function $(r, ^, z) in Eq. (2) into 

<^{r,<P,z) = R,{r)e^""^Z{z) 
leads to the longitudinal wave equation 

J2 



dz^-' 



-\- m 



*2 



Z{z] 







(3) 



and the radial wave equation 




dr"^ r dr 



RJr) 



(4) 



where C.^ 



m 



m^(r), and m* is separation constant which physically is the transverse 



mass. For simplicity we shall study the case of an infinitely long cylinder; this gives E'^ = 
P^__ -\- m*^, where E is the eigenenergy and {PnJ\ S'l'e the eigenvalues associated with the 
longitudinal wave functions Zn^{z) = Ae^'^"^^ . Eq. (4) is a Bessel equation which has the 
following two solutions: 

Rl{r) = AJiy{Ciyr) in region I, 



R 



III 



A'K^[C'lr) in region II . 



where C.^ 



m 



*2 



m^, C. 
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m 



*2 



and A and A' are overall normalizations. Using 



the continuity of the logarithmic derivatives of the two wave functions Rl{r) and R 



III 



and taking the asymptotic forms of both Jj, and Ky^ as M 
eigenvalue equation for Ci, as 



Ru{Ro) — Jiy{Ciy^sRo 







00 in region II, we obtain an 



(5) 
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which gives C,y^s = iv,sl Ro^ where the values {jv^s} are the zeros of the Bessel function Jj,. 
The radial wave function is then 

We note that for the case of transverse confinement, m* is quantized and contributes a 
zero-point energy associated with motion in the transverse direction. Thus, as is evident 
from Eq. (3), even for a gluon with zero rest-mass, the gluon acquires a non-zero intrinsic 
momentum distribution due to the kinetic energy of confinement. 

Introducing the Fourier transforms of Zn^{z) and R^^s find longitudinal and 

transverse wave functions in momentum space, 

= 4= / e^'^^- Z^^{z)dz (6) 
v Ztt J 



anc 



2tt 

5_L I rdr Mp^r)MC,,sr) , (7) 

ZTT 



where i? is a constant. Carrying out the integration gives [18] 

R 

^J'.sIPt) = 7^ [Pt Ji,+i{PtRo) JACi,,sRo] 

~Ciy,s Ju{PtRo) Ju+i{Cu,sRo)\ 



where 



R 

\Pt Ju+i{PtRo) Ju{Cu,sRo 



2 _ r<2 

Pt ^l': 

~Ciy,s Ju{PrpRo) Ju+l{Cu,sRo)\ 



The wave function ^j/,s(Pt) varies approximately as p^"^, with a small oscillatory component, 
at large . The single-particle wave function in momentum space is defined as 

and hence the gluon momentum distribution (Eq. (1)) becomes 

where the gluon degeneracy number (/^ is 8. The transverse momentum distribution of the 
gluons in a quark-gluon plasma of temperature T and chemical potential // is therefore given 
by, using Eqs. (6), (8), and (9), 
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J3p(Jz/27r) 



R 



p2 - C2 



■ (10) 



In contrast, if the boundary of the quark-gluon plasma extends to infinity, then the 
transverse momentum distribution is given by the Bose- Einstein form: 



J3p(Jz/27r) (27r)2 eiE-^)IT _ i 

Fig. 1 shows the transverse momentum distributions of the gluons in a quark-gluon 
plasma of different radii Ro = 2, 3.5, and 5 fm, at temperature T = 200 MeV, // = 0, and 
Pz = The solid curves give the gluon transverse momentum distributions including the 
boundary effects, as given by Eq. (10). The dashed curves give those distributions for an 
unconfined plasma, obtained using the Bose- Einstein distribution, given by Eq. (If). Eor a 
plasma of a small radius, the distribution given by Eq. (10) (the solid curves) differs from 
that given by the Bose- Einstein distribution of Eq. (II) (the dashed curves), as shown in 
Eig. 1(a). In each case, the distribution which includes the boundary effect is suppressed at 
low relative to the Bose- Einstein distribution; this is due to the small density of confined 
single-particle states at low energies. As the radius of the plasma increases, the distribution 
in the low momentum region increases and approaches the Bose- Einstein distribution of an 
unbounded boson system; this is shown in Eig. 1(b) and Eig. 1(c). The large difference 
between the two distributions at very low momenta, p^ < O.I GeV, is due to the fact that 
Bose- Einstein distribution diverges as the momentum of the massless gluon goes to zero, and 
in contrast the distribution with a boundary does not, due to the smaller \ow-p^ density of 
states. In the opposite limit of large p^, the bounded case has a large high-momentum tail, 
which varies with the transverse momentum by a power law, and exceeds the Bose- Einstein 
distribution by many orders of magnitude. This is due to the vanishing of the single-particle 
wave functions at the boundary of the plasma. The momentum at which the distribution 
begins to differ substantially from the Bose- Einstein distribution increases with the radius 
of the plasma. 

In contrast to our results for the momentum distribution at low-p^, Sinyukov [19] studied 
the finite-size effect on the momentum distribution of a massless particle using Bjorkan's 
scale hydrodynamics. He found that the single-particle transverse momentum at p^, 
is much enhanced over the Bose- Einstein distribution (The B.-E. distribution is infinite at 
for an equilibrated massless particle gas). The physical reason of the enhancement 
is not known. His results also contradicts those of Bolz et al. [17] and Kataja et al. [20] 
who studied the transverse momentum distribution using Bjorkan's scale hydrodynamics. 
The p^ distributions for massless particles in the latter two calculations do not show an 
enhancement over the B.-E. distribution. 



IV. EFFECT OF MEAN-FIELS BOUNDARY ON A PION SYSTEM 

After an ultrarelativistic nucleus-nucleus collision the produced pions are at high tem- 
perature, typically 100 MeV to 200 MeV. It is well known that tt-tt interaction in the isospin 
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/ = state is attractive [21]. Furthermore, the strong interaction range of the pions is about 
1.4 fm and the average separation between pions, J, in these pion systems can be estimated 
as follows. The pion density is approximately [22,5] 

1 dN^ 
n = — — — 
A tody 

where to = I fm/c is the pion production time and A is the transverse area of the collision 
region. For ^^0 on Au at 200A GeV, the peak value of dN^/dy is 120 • | = 180 and it 
is 120 for S on S at 200A GeV. The average separation, d = between pions in the 

pion system produced in the former collision is about 0.6 fm and about 0.7 fm in the latter 
collision. Because the average separation between pions is smaller than the pion interaction 
range, the pions in such systems are strongly interacting with each other and they are moving 
in an attractive mean-field potential which extends over the pion system. The attractive 
mean-field leads to a quasi-bound pion system which can be considered to be in a liquid 
phase with a surface boundary, rather than in a gas phase. The notion of a pion liquid was 
introduced by Shuryak [23], to emphasize the importance of interparticle interactions; these 
interactions are attractive and have easily observable consequences. 

We assume that hadronic matter after a heavy-ion collision is only confined within the 
boundary of their mean field when their average separation is smaller than the interaction 
range between pions. They become free pions when the pions' average separation in the 
system become larger than their interaction range. Because of the short range nature of 
the pion interactions, as the system expands, the transition from a pion with an attractive 
mean field potential to a free pion system, i.e. freeze-out, is a rapid process. The momentum 
distribution of the pions is governed much by their momentum distribution just before they 
freeze-out. It is therefore reasonable to use the boundary of their mean field before the freeze- 
out transition to calculate the momentum distribution of the observed pions. Pions in such 
a system have a momentum distribution that is modified by the presence of this boundary. 
For example, the boundary leads to a vanishing of the pion wave function throughout the 
exterior region, which enhances the momentum distribution in the high momenta region 
relative to a Bose-Finstein distribution. It is interesting to apply the previous formulation 
to this pion system, to develop our understanding of the boundary effects on a massive boson 
system. 

We consider that a pion system in a cylinder with a sharp transverse (radial) boundary. 
We envision the boundary disappearing at the moment of freeze-out, while the momentum 
distribution of pions remain unchanged. The mean-field potential m(r) for this pion system 
(before freeze-out) is taken to be 

, . I m if r < Ro (region I) 

Mr) = { " -r ~ D / TTX 

I oo n r > Ko (region iij , 

where is the physical pion mass. 

Following the discussion of section III, the pion transverse momentum distribution in 
the bounded system, with temperature T and chemical potential is given by 

J3p(J^/2vr) ^ L J 
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R 



2 _ r<2 

Pt ^i'. 



12) 



In contrast, the pion transverse momentum distribution without a boundary is given by 
the Bose- Einstein distribution, 

dN„ ttRI 1 

J3p(Jz/27r) (27r)2 e(^-^)/^ - 1 ' ^ ' 

where E"^ = + + and = 3. 

Fig. 2 shows the pion transverse momentum distributions arising with the boundary 
(using Eq. (12), sohd curves) and with a Bose- Einstein Distribution (Eq. (13), dashed curves) 
for systems of different radii: Ro = 2,4, and 6 fm, T = 200 MeV, Pz = 0, and // = 0. The 
general features of the pion transverse momentum distribution, shown in Eig. 2, are similar 
to that of the gluons, shown in Eig. 1. Eor a pion system of small radius, the distribution 
arising from the boundary effect (solid curves) is quite different from that of a Bose- Einstein 
distribution (dashed curves) but they approach each other as the size of the system goes to 
infinity. Note however that the pion transverse momentum distribution at low momenta is 
broader than we found for gluons. 

Eig. 3 shows the pion transverse momentum distribution Eq. (12) for a pion system 
with a boundary (solid curve), and the Bose- Einstein distribution Eq. (13) (dashed curve), 
compared with experimental data for S-|-S — > negatives -\- \ at E = 200A GeV [6-8]. Eig. 4 
similarly shows 0-|-Au — )• negatives -\- X . The radius Ro of the pion source is taken from 
the results of pion interferometry [24]. By equating the root-mean-square radius of our sharp 
distribution and the Gaussian distribution of [24], we have Ro = \/2R^^ where R^ is the 
effective transverse size of the source, as defined in [24]. Our calculations fit the experimental 
data well over most of p^, given a pion system of radius Ro = 5.0 fm for S-|-S and Ro = 7.0 
fm for 0-|-Au, at T = 200 MeV, // = 0, and = 0. However, at low momenta [p^ < m^) 
neither the bounded pion system distribution Eq. (12) (solid curve) nor the Bose- Einstein 
distribution Eq. (13) (dashed curve) fits the data. 



V. COULOMB EFFECTS ON A PION SYSTEM 

The Coulomb final-state interaction in relativistic nuclear collisions has been studied by 
many authors [25-27]. This effect has been found to play an important role in the momen- 
tum distribution of charged multiplicity produced in relativistic nucleus-nucleus collisions. 
Eor example, Libbrecht and Koonin [26] used a covariant classical formulation to show that 
a Coulomb focusing could arise toward Oc.m. = 90° and finite p^, in the Coulomb field gener- 
ated by two charged fragments in relative motion (the projectile and the target remnants). 
Gyulassy and Kauffmann [27] used a relativistic field theoretic model to derive formulas to 
first order in Za for the Coulomb final-state interactions. They derived covariant formu- 
las that took into account multiple independently moving charged fragments of finite size 
and finite thermal expansion velocities, and showed the importance and complexity of the 
Coulomb final-state interaction in nuclear collisions. 
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The charged pion inclusive cross-section (T^±(p), where p is the particle momentum 
to first order in Za and to first order in the gradient of the pion-source current J[x) is given 
(by Eq. (3.13) in [27]) as 

a^±{p) ^ aAp T Sk{p)){l T 8D{p)] , (14) 

where cr^o(p) is the neutral pion inclusive cross-section, 6k[p) is the shift in the momentum^ 
due to the Coulomb final-state interaction, {1 =F ^D[p)} is the Coulomb phase space factor, 
and the it signs refers to the charge of the pion considered. In our calculations cr^o(p) is 
identified either by the distribution of a bounded pion system, Eq. (12), or by a Bose- Einstein 
distribution, Eq. (13). The longitudinal momentum p^, in terms of the rapidity y and the 
transverse mass m^,, is 

p_, = sinh y , (15) 



-|- . (16) 

The Coulomb phase space factor for a static charge density and in the limit of small mo- 
mentum is found to be given by the Gamow factor G{r])^ which is defined as: 

where for a ib|e| charged particle, rj = ±Za/v^ where Z is the total charge, a is the fine 
structure constant, and v is the pion velocity. In our calculations we replace the first order 
expansion {1 =p 6D} of Eq. (14)by the Gamow factor G{r])^ to include Coulomb effect to all 
orders of Za. Eq. (14) thus becomes 

a^±ip)^aApTSkip))-^^^. (18) 

Eig. 5 shows how the Coulomb final-state interaction given by Eq. (18) affects the pion 
transverse momentum distribution, which was obtained previously using Eq. (12) for a pion 
system formed after a relativistic S-|-S collision, with Z = 32 and Ro = 5.0 fm. It is clear 
that, at low p^, the Coulomb final-state interaction enhances the 7r~ transverse momentum 
distribution and suppresses that of tt"'" relative to tt". 

Eig. 6 shows, relative to that of a p -\- p collision in comparable energies, the ratios of 
the normalized pion transverse momentum distributions for a bounded pion system given 
by Eq.(12) (solid curve), and the Bose- Einstein distribution of Eq.(13) (dashed curve), com- 
pared to experiment for S-|-S — > negatives -\- X [28]. The Coulomb final-state interaction 
defined by Eq. (18) is taken into account, and calculations are carried out in the rapidity 
range 2.0 < y^^^ < 3.0 for a system of radius Ro = 5.0 fm, with Z = 32, T = 200 MeV, 
/i = 0, and pz = 0. The transverse momentum distribution of the charged pions is found to 
be quite sensitive to the rapidity of the pion. The pion transverse momentum distribution 
from p -\- p collisions is parametrized by [28] 

1 1 dN J_ /„2 1^2 

— ^ e-BVfT+^ , 19 
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where B = 0.162 it 0.003 GeV and C = 0.11 it 0.02 GeV, and the normahzation is defined 
by 

.2 Gey dN ,^ 
/ -J^dp^ = 1- 

Fig. 7 shows the ratios of the normahzed pion transverse momentum distribution relative 
to that from p -\- p coUisions at comparable energies for a pion system with a boundary 
(solid curve) and for the Bose- Einstein distribution (dashed curve), in comparison with 
experimental data for 0+Au — )• negatives + X [28]. The Coulomb final-state interaction 
given by Eq. (18) is again included, and calculations were carried out for the rapidity range 
2.0 < y^^^ < 3.0 for a system of radius Ro = 7.0 fm, with Z = 60, T = 200 MeV, // = 0, and 

Pz = 0/ 

The last two figures, Eig. 8 and Eig. 9, show a direct comparison of these pion transverse 
momentum distributions, for the case with a boundary (Eq. (12), Eigs. 9 (a) and 8 (a)) 
and Bose- Einstein (Eq. (13), Eigs. 9 (b) and 8 (b)) with data from S-|-S — > negatives -\- X 
and 0-|-Au — )• negatives -\- X [6-8], taking into account the Coulomb final-state interaction. 
Both sets of data were fitted using T = 200 MeV, // = 0, 2.0 < y^^^ < 3.0, and Ro = 5.0 
fm and Z = 32 for S-|-S data and Ro = 7.0 fm and Z = 60 for 0-|-Au data. It is clear 
that incorporation of both the boundary effects and the Coulomb final-state interaction 
effects on the pion system leads to good agreement with the data over the range of p^, under 
consideration, as shown in Eigs. 9 (a) and 8 (a). Also, Bose- Einstein distribution along with 
Coulomb final-state interactions give good agreement with data and it give a slightly larger 
enhancement at low-p^, as in Eigs. 9 (b) and 8 (b). 

Coulomb final-state interactions predict also that there are significant differences between 
the transverse momentum distributions of tt"'", tt", and tt" in the low transverse momentum 
region. 



VI. SUMMARY AND CONCLUSIONS 

The effect of a boundary on a system of massless bosons (gluons) and a system of massive 
bosons (pions) at high temperature have been studied. The boundary is found to affect the 
transverse momentum distribution of these systems. Relative to a free-particle Bose- Einstein 
distribution, the transverse momentum distribution of gluons and pions is suppressed at low 
momenta and enhanced at high momenta. Eor the gluons in a quark-gluon plasma of a 
radius of only a few times T"^, the transverse momentum distribution at low momenta 
relative to R~^ is much smaller than the Bose- Einstein distribution and approaches it as the 
radius of the plasma goes to infinity. This suppression is due to the lower level density of 
single-particle states at the bottom of the well for a finite system [29]. It can also be depicted 
as due to the weakness of large- wavelength amplitudes in a bound system. In contrast, the 
high momentum tail arises from the rapid variation of the single-particle wave function at 
the boundary of the system. It is clear that the boundary effects on a quark-gluon plasma 
play an important role in the gluon transverse momentum distribution, and consequently 
on signatures of the plasma. The boundary effects should be taken into account to provide 
more reliable estimates of the signatures of quark-gluon plasma-formation. 

The pions in a hadronic system produced after a high energy nucleus-nucleus collision 
interact with each other and with the surrounding hadronic and nuclear matter. These 
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interactions, in the isospin state / = 0, are attractive. The pions can be modeled as moving 
in a mean-field potential which creates an effective boundary for the system. Here we 
assume that the pion system has a sharp boundary before freeze-out. The general features 
of the pion transverse momentum distribution are rather similar to that of the gluons. 
There is a suppression at low momenta, with a scale that depends on the radius of the 
system, and there is a large enhancement relative to Bose- Einstein distribution at high 
momenta, p^, > 2.0 GeV. The distribution is quite different from Bose- Einstein for systems 
of small radius, and approaches Bose- Einstein as the radius is taken to infinity. However, 
the pion transverse momentum distribution is broader than that of gluons at low momenta. 
Eor > m^, good agreement with experimental data is obtained for the pion transverse 
momentum distributions for S-|-S — > negatives -\- X and 0-|-Au — )• negatives -\- X . 

The pion transverse momentum distribution at low p^, for a system with a boundary does 
not depart considerably from a Bose- Einstein distribution for large systems, when compared 
with the available experimental data. However, the bounded model predicts that there is 
always an enhancement relative to the Bose- Einstein distribution at high momenta, p^, > 2.0 
GeV. There is a possible contamination from hard scattering processes to the high transverse 
momentum tail, which makes the boundary effect difficult to observe. However, structure 
function is poorly known in this p^, region and a quantitative study is needed in the future 
to differentiate between the two mechanisms. 

What is the momentum distribution of an expanding pion system. The transverse mo- 
mentum distribution of the pions is suppressed at low p^, (relative to static system) giving 
rise to a depression in the distribution in that region, if the particles in the system have been 
boosted by the same amount. However, if the boost increases with the energy of the single 
particle wave function, a very small suppression in the transverse momentum distributions 
at low is found and the magnitude of the suppression is slightly increases with the amount 
of the boost. On the other hand, we have found no change in the high transverse momen- 
tum tail with the magnitude of the boost, while the intermediate portion of the distribution 
slightly increases. A pion system with a boost behaves as if it is approximately a system 
with a higher temperature. 

The puzzle of the pion transverse momentum distribution in the low p^, region, p^, < m^, 
where an enhancement is observed in the data, leads us to study the effect of Coulomb 
final-state interactions on the system. Coulomb final-state interaction on pions arise from 
target and projectile remnants near the central rapidity region. This Coulomb final-state 
interaction modifies the charged pion transverse momentum distributions considerably at 
low transverse momenta, where the distribution of charged pions relative to neutral pions is 
enhanced for 7r~ but suppressed for tt"'". The degree of enhancement or suppression is quite 
sensitive to the rapidity of the particles in the system. The inclusion of effects calculated 
from Coulomb final-state interactions leads to very good agreement with 7r~ experiment, 
differences between the transverse momentum distributions of tt"'", tt", and tt" in the low 
transverse momentum region await future experimental data for direct comparison. 
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FIGURES 



FIG. 1. The transverse momentum distributions of gluons in a quark-gluon plasma with a 
sharp cylindrical boundary, for radii Rg = 2, 3.5, and 5 fm, at temperature T = 200 MeV, = 0, 

= 0, and = 0. The solid curves are the numerical results obtained using Eq. (10), and the 
dashed curves are the predictions for a Bose- Einstein distribution, given by Eq. (11). 

FIG. 2. The transverse momentum distributions of pions in a pion system with a sharp cylin- 
drical boundary, for radii Rg = 2, 4, and 6 fm, at temperature T = 200 MeV, = 0, = 140 
MeV, and = 0. The solid curves are the distributions arising from the boundary effect, given by 
Eq. (12). The dashed curves are the predictions for a Bose- Einstein distribution, given by Eq. (13). 

FIG. 3. The pion transverse momentum distributions in a pion system with a boundary given 
by Eq. (12) (solid curves), and the predictions from the Bose- Einstein distribution given by Eq. (13) 
(dashed curves), in comparison with the experimental data for S + S negatives -|- X at E = 200A 
GeV [6-8]. Good agreement is obtained with a radius Rg = 5 fm, temperature T = 200 MeV, 
^ = 0, = 140 MeV, and = 0. 

FIG. 4. The pion transverse momentum distributions from Eq. (12) for a bounded pion system 
(solid curves), and a Bose- Einstein distribution (Eq. (13), dashed curves), in comparison with 
experiment for 0-|-Au negatives -|- X at E = 200 A GeV [6-8]. Good agreement is obtained 
with a radius Rg = 7 fm, temperature T = 200 MeV, = 0, = 140 MeV, and p^ = 0. 

FIG. 5. The effect of the Coulomb final-state interaction, calculated using Eq. (18) on the 
transverse momentum distributions of vr"*", vr", and7r° in a bounded pion system produced after the 
collision S + S negatives + X (see text). 

FIG. 6. The ratios, relative to that of a collisions, of the normalized transverse momentum 
distributions of negatively charged pions in a pion system with a cylindrical boundary of radius 
Rg = 5.0 fm, given by Eq. (12) (solid curve), and the Bose- Einstein distribution, Eq. (13) (dashed 
curve), in comparison with data from the collision S + S negatives + X [28]. For the bounded 
pion system results have been calculated taking into account Coulomb final-state interaction, as in 
Eq. (18) (see text). 

FIG. 7. The ratios of the normalized transverse momentum distributions of negatively charged 
pions in a pion system with a cylindrical boundary of radius Rg = 7.0 fm (Eq. (12), solid curve), 
and the Bose- Einstein distribution (Eq. (13), dashed curve), relative to that of p -|- p collisions, 
in comparison with data from the collision 0-|-Au negatives + X [28]. For the bounded pion 
system results have been calculated taking into account Coulomb final-state interaction, as in 
Eq. (18) (see text). 
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FIG. 8. The pion transverse momentum distributions in a pion system with a boundary, 
Eq. (12), as in Eq. (18) (fig. (a)), and the predictions from Bose- Einstein distribution, Eq. (13) 
(fig.(b)), taking into account Coulomb final-state interaction, in comparison with experiment for 
S + S negatives + X [6-8]. Theoretical results are obtained with radius Rg = 5 fm, for Z = 32, 
temperature T = 200 MeV, = 0, = 140 MeV, and 2.0 < y < 3.0. 

FIG. 9. The pion transverse momentum distributions in a pion system with a boundary, 
Eq. (12), as in Eq. (18) (fig. (a)), and the predictions from Bose- Einstein distribution, Eq. (13) 
( fig.(b)), taking into account Coulomb final-state interaction, in comparison with experiment 
for 0-|-Au negatives + X [6-8]. Theoretical results are obtained with radius Rg = 7 fm, for 
Z = 60, temperature T = 200 MeV, = 0, = 140 MeV, and 2.0 < y < 3.0. 
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